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ABSTRACT 



Aims. To determine alpha effect and turbulent magnetic diffusivity for mean magnetic fields with profiles of different length scale 
from simulations of isotropic turbulence, and to relate these results to nonlocal formulations in which alpha and the turbulent magnetic 
diffusivity correspond to integral kernels. 

Methods. A set of evolution equations for magnetic fields is solved which gives the response to imposed test fields, that is, mean 
magnetic fields with various wavenumbers. Both an imposed fully helical steady flow consisting of a pattern of screw-like mo- 
tions (Roberts flow) and time-dependent statistically steady isotropic turbulence are considered. In the latter case the aforementioned 
evolution equations are solved simultaneously with the momentum and continuity equations. The corresponding results for the elec- 
tromotive force are used to calculate alpha and magnetic diffusivity tensors. 

Results. For both the Roberts flow under the second-order correlation approximation and isotropic turbulence alpha and turbulent 
magnetic diffusivity are largest at large scales and these values diminish toward smaller scales. In both cases the alpha effect and tur- 
bulent diffusion kernels are well approximated by exponentials, corresponding to Lorentzian profiles in Fourier space. For isotropic 
turbulence the turbulent diffusion kernel is half as wide as the alpha effect kernel. For the Roberts flow beyond the second-order 
correlation approximation the turbulent diffusion kernel becomes negative at large scales. 

Key words. Magnetohydrodynamics (MHD) - Turbulence 
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1. Introduction 

Stars and galaxies harbour magnetic fields whose scales are large 
compared with the scale of the underlying turbulence. This phe- 
nomenon is successfully explained in terms of mean-field dy- 
namo theory discussed in detail in a number of text books and re- 
views (e.g. Moffatt 1978, Krause & Radler 1980, Brandenburg & 
Subramanian 2005a). In this context velocity and magnetic fields 
are split into large-scale and small-scale components, U -U+u 
and B = B + b, respectively. The crucial quantity of the the- 
ory is the mean electromotive force caused by the small-scale 
fields, £ = H X ft. In many representations it is discussed un- 
der strongly simplifying assumptions. Often the relationship be- 
tween the mean electromotive force and the mean magnetic field 
is tacitly taken as (almost) local and as instantaneous, that is, £ 
in a given point in space and time is considered as determined 
by B and its first spatial derivatives in this point only, and the 
possibility of a small-scale dynamo is ignored. Then the mean 
electromotive force is given by 



TjijkdBjIdxk 



(1) 



with two tensors and riij]^. If the turbulence is isotropic the 
two tensors are isotropic, too, that is aij - adij and rjijk = rjteijic 
with two scalar coefficients a and rj^. Then the expression ([T]l 
simplifies to 



£ = aB - T]iJ , 



(2) 



where we have denoted V x B simply by / (so that / is yUo times 
the electric current density, where yUo is the magnetic permeabil- 
ity of free space). The coefficient a is, unlike //t, only non-zero 



if the turbulence lacks mirror-symmetry. The coefficient rji is re- 
ferred to as the turbulent magnetic diffusivity. 

In general, the mean electromotive force has the form 



G^Go + KoB, 



(3) 



where £o stands for a part of £ that is independent of B, and 
K o B denotes a convolution in space and time of a kernel K 
with B (see, e.g., Krause & Radler 1980, Radler 2000, Radler & 
Rheinhardt 2007). Due to this convolution, £ in a given point in 
space and time depends on B in a certain neighborhood of this 
point, with the exception of future times. This corresponds to a 
modification of ([T]i such that also higher spatial and also time 
derivatives of B occur. 

In this paper we ignore the possibility of coherent effects 
resulting from small-scale dynamo action and put therefore £o 
equal to zero. For the sake of simplicity we assume further the 
connection between £ and B to be instantaneous so that the con- 
volution K o B refers to space coordinates only. The memory 
effect, which we so ignore, has been studied previously by solv- 
ing an evolution equation for £ (Blackman & Field 2002). 

For homogeneous isotropic turbulence we may then write, 
analogously to Q, 



£>-aoB-%oJ, 

or, in more explicit form. 



£(^) = J [a(mx UOJix - ^)] d'^ 



(4) 



(5) 
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with two functions a and ?7t of ^ = |^| that vanish for large ^. The 
integration is in general over all ^-space. Although £ and B as 
well as a and fji may depend on time the argument t is dropped 
everywhere. For a detailed justification of the relations (|4| and 
(|5]) we refer to AppendixlAl In the limit of a weak dependence of 
B and J on space coordinates, i.e. when the variations of B{x—^) 
and Jix - f ) with ^ are small in the range of ^ where a{^) and 
77t(^) are markedly different from zero, the relations (|4]l or (|5]l 
turn into (|2]i, and we see that a - ^ a{^) d?^ and Jjt - J fjii^) d^*^. 

At first glance the representations (|4]l and (|5]) of £ look rather 
different from (|3]l. Considering J - V x B and carrying out an 
integration by parts we may however easily rewrite (|5]l into 



£,W = J Kij(^Bj(x-^)d'^ 
with 



(6) 



(7) 



We further note that due to the symmetry of a{^) in ^ only the 
part of B{x - ^) that is symmetric in ^, i.e. the part that can be 
described by B{x) and its derivatives of even order, contributes 
to the a terms in (|5]) or in ^ and (|7]l. The symmetry of ffti^) 
implies that only the part of B{x - ^) antisymmetric in ^, which 
corresponds to the derivatives of B(x) of odd order, contributes 
to the f)t terms. 

Finally, referring to a Cartesian coordinate system (x,y,z) 
we define mean fields by averaging over all x and y so that in 
particular £ and B depend only on z and on time. Then (|5]l turns 
into 



8(z) = J [a(OB(z -0- UOJiz - O] . 



(8) 



The functions a(0 and %(0 are just averages of a(^) and fji(^) 
over all and ^y. They are therefore real and symmetric in = 
^. The integration in (O is in general over all The remark on 
the limit of weak dependences of B and J on space coordinates 
made above in connection with (|4|l and (|5]l applies analogously 
to (|8]l. We have now a - J a(l^) d^ and rjt - J fltiO d(- 

Relation ([8]l can also be brought in a form analogous to @ 
and (I2]l, 



-Jk, 



with 



QUO 



(9) 



(10) 



The remarks made under (|6]l and (|7]) apply, now due to the sym- 
metries of a{^) and ^t(^) in analogously to ([8]l, (|9]l and ( fTOl l. 

It is useful to consider in addition to ([8]l also the correspond- 
ing Fourier representation. We define the Fourier transformation 
in this paper by Q{z) — J Q(k) exp(ife) d(k/2K). Then this repre- 
sentation reads 



G{k) = aik)Bik) - Uk)J(k) . 



(11) 



Both a(k) and ^t(^) are real quantities, and they are symmetric in 
k. The limit of weak dependences of B and J on z corresponds 
here to A; — > 0, and we have a - a(Q) and T/t = ^t(O). Detailed 
analytic expressions for a(f) and fjtiO, or aik) and ^t(^), can be 



derived, e.g., from results presented in Krause & Radler (1980). 
A numerical determination of quantities corresponding to a(0 
and fjtiO has been attempted by Brandenburg & Sokoloff (2002) 
for shear flow turbulence. 

In this paper two specifications of the velocity field u will be 
considered. In the first case u is chosen such that it corresponds 
to a steady Roberts flow, which is periodic in x and y and inde- 
pendent of z- A mean-field theory of a magnetic field in fluid 
flows of this type, that are of course different from genuine tur- 
bulence, has been developed in the context of the Karlsruhe dy- 
namo experiment (Radler et al. 2002a,b, Radler & Brandenburg 
2003). It turned out that the mean electromotive force £, except 
its z component, satisfies relation ^ if any nonlocality in the 
above sense is ignored (see also AppendixlBl). Several analytical 
and numerical results are available for comparison with those of 
this paper In the second case u is understood as homogeneous, 
isotropic, statistically steady turbulence, for which the above ex- 
planations apply immediately. Employing the method developed 
by Schrinner et al. (2005, 2007) we will in both cases numeri- 
cally calculate the functions a{k) and f}i{k) as well as a{() and 
^t(0- 

2. The method 

2.1. A more general view on the problem 

We first relax the assumption of isotropic turbulence used in the 
Sect.[T](but will later return to it). We remain however with the 
definition of mean fields by averaging over all x and y. Then, as 
already roughly indicated above, and By may only depend on 
z and time but B,, because of V ■ B = 0, must be independent 
of z. Furthermore all first-order spatial derivatives of B can be 
expressed by the components of V x B, that is, of J, where = 
0. Instead of ^ we have then 



£iiz) = j [&ij(OBj(z -0- flijiOJjiz - O] 

and instead of (fTTl l 

%{k) = aij(k)%{k) - iiij(k)lj(k) , 



(12) 



(13) 



with real and fjiji^), which are even in and real atjik) 

and fjijik), which are even in k. A justification of these relations 
is given in AppendixlAl We have further 

&ij{k)^ j aijiO cos k(d(, fjij(k) = j fiijiO cos k(d(. (14) 

Since 73 = the fjn, as well as the 77,3, are of no interest. 

In the following we restrict attention to £f and £,. and assume 
that B, is equal to zero. We note that £, and the contributions of 
to £v and £,. are anyway without interest for the mean-field 
induction equation, which contains £ only in the form V x £, 
that is, they do not affect the evolution of B. We may formulate 
the above restriction in a slightly different way by saying that we 

consider in the following £,, aij and rjij as well as £,, aij and fjij 
only for 1 < j < 2. 

As for the mean-field treatment of a Roberts flow depend- 
ing on X and y only (and not on z) we refer to the aforemen- 
tioned studies (Radler et al. 2002a,b, Radler & Brandenburg 
2003). Following the ideas explained there we may conclude that 
- &i0^ij and fjiji^) - fjtiO^ij with functions a and fjt of 
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^, and analogously aij(k) = a(k)6ij and fjijik) - fit{k)6ij with 
functions a and fjt of k, all for 1 < j < 2. For obvious reasons 
the same is true for homogeneous isotropic turbulence. 



From ( fT9] l and (|20] | we conclude 



da — — 

— -Uxb+uxB+uxb 

at 



u X b - rjj , 



(21) 



2.2. Test field method 

We calculate the aij(k) and i)ij(k), or a{k) and fi{k), numerically 
by employing the test field method of Schrinner et al. (2005, 
2007). It has been originally developed for the calculation of the 
full a and tj tensors [in the sense of ([T]l] for convection in a spher- 
ical shell. Brandenburg (2005) employed this method to obtain 
results for stratified shear flow turbulence in a local cartesian do- 
main using the shearing sheet approximation. More recently, Sur 
et al. (2008) calculated in this way the dependences of a and rfi 
for isotropic turbulence on the magnetic Reynolds number, and 
Brandenburg et al. (2008) have calculated the magnetic diffusiv- 
ity tensor for rotating and shear flow turbulence. However, in all 
these cases no nonlocality in the connection between £ and B 
has been taken into account. 

Following the idea of Schrinner et al. we first derive expres- 
sions for £ with several specific B, which we call "test fields". 
We denote the latter by B and definqj 



b'' = B(cosfe,0,0), 



b'' 



B (0, cos kz, 0), 
b'' = B(smkz,Q,Q), B^' = B(0,sinfe,0), 



(15) 



with any constant B and any fixed value of k. We then replace 
B and 7 in (O by b'"' and V x b' "^ or by b'" and V x b''\ 

p C p S 

Denoting the corresponding Shy S or by S , respectively, 
and using (fT4l > we find 



£ (z) 



B ip(k) cos kz - ri\^ik)k sin kz^ , 
£f (z) - B ^aip(k)smkz + ii]^j(k)kcoskz\ , 
for 1 < i, p <2, where 

From this we conclude 



(16) 



(17) 



an(k) = B 



-1 



— 7 c — 7 s 

£; (z) COS kz + &j (z) sin kz 



f,\.{k) = -{kBY 



£■ '^(z) sin kz - &\^{z) cos kz 



(18) 



for 1 < /, < 2. 

These relations allow us to calculate the Uij and rjij if the 

— 

&i with 1 < i,p < 2 for both q - c and ^ = s are known. In 
preparing the numerical calculation we start from the induction 
equation. Its uncurled form reads 



'4-UxB. 

dt 



rjJ, 



(19) 



where A is the magnetic vector potential, B = V x A, and J - 
V X B. Here the Weyl gauge of A is used. Taking the average of 
( fT9] l we obtain 



QA. _ _ _ 

— - UxB + uxb-JiJ. 

ot 



(20) 



' The notation used here differs slightly from that in Brandenburg et 
al. (2008), where first test fields B were introduced and only later the 
two versions B and B are considered. 



where a - A - A and j - J - J -V xb. 

For the calculation of the e1 ^ we are interested in the 6'"' = 
V X a'"' which occur in response to the test fields B . Specifying 
(l2n i in that sense we obtairQ 



daPi - 



dt 



^UxbP'' +uxB" +uxb'"' -uxb'"' - Tjf 



■pq 



(22) 



Equations of this type are called "test field equations". 

So far no approximation such as the second order correla- 
tion approximation (SOCA), also known as first order smooth- 
ing approximation, has been made. If we were to make this as- 
sumption, terms that are nonlinear in the fluctuations would be 
neglected and ( l22l i would simplify to 



^^Uxb"^+uxB'" 
dt 



rjf' (for SOCA only). 



(23) 



In the following SOCA results will be shown in some cases for 
comparison only. 

In either of the two cases the aij and fiij are to be calculated 



from £'* = M X b'"' . More details of the numerical calculations 



of the £ will be given below in Sect. 12.41 

Returning once more to ( 1181 ) we note that the £ depend 

— PQ 

on both k and z introduced with the B . As a consequence of 
imperfect simulations of the turbulence they may also depend 
on the time t. The and rjij however should depend on k but 
no longer on z and f. We remove the latter dependences of our 
results by averaging S/j and rjij over z and t. For the Roberts flow 
there should be no such z or f dependences. 

The relations ( fTSl l allow the determination of all components 
of Qij and Tjij with 1 < /, y < 2. We know already that aij - adij 
and fjij = fjtSij, that is, an = 522 = a, ^ii = ?722 = ^ and 
a\2 - Q'21 - fjii = Jjii - 0. We may therefore determine a and 
according to a = an and iit - fjn by using the two test fields 

B * and the relations (fTSl l with / = j - ^ only. 

2.3. Flow fields 
2.3.1. Roberts flow 

We consider here a special form of a steady flow which, in 
view of its dynamo action, has already been studied by Roberts 
(1972). It has no mean part, [7 = 0, and u is given by 



u - -zxVip + kfipz , 
where 

— (uQ/ko)coskox coskoy, kf — y/2ko 



(24) 



(25) 



with some constant ko. The flow is fully helical, V x m = k[U. 
The component form of u as defined by jTM and (|25] | reads 

u = Mo(-cosA;ox sin k^y, sin kox coskoy, V2cosA:ox cosA:o3').(26) 
We note that = u^. 



^ In the corresponding expression (27) of Brandenburg (2005) the U 
term is incorrect. This did not affect his results because U = 0. 
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2.3.2. Turbulence 

Next, we consider isotropic, weakly compressible turbulence 
and use an isothermal equation of state with constant speed of 
sound, Cs- Considering first the full velocity field U - U + u we 
thus accept the momentum equation in the form 



dU 
dt 



— = -I/- V[/-c^Vlnp+/+p"'V-2pvS, 



(27) 



where / is a random forcing function consisting of circularly po- 
larized plane waves with positive helicity and random direction, 
and S/j = j( Uij+Ujj)--^6ij'V-U is the traceless rate of strain ten- 
sor The forcing function is chosen such that the magnitude of the 
wavevectors, |, is in a narrow interval around an average value, 
which is simply denoted by kf. The corresponding scale, k^^, is 
also referred to as the outer scale or the energy-carrying scale of 
the turbulence. More details concerning the forcing function are 
given in the appendix of Brandenburg & Subramanian (2005b). 
With the intention to study the mean electromotive force in the 
purely kinematic limit the Lorentz force has been ignored. 

In addition to the momentum equation we use the continuity 
equation in the form 



dlnp 
dt 



= -U-Vlnp-V-U. 



(28) 



In all simulations presented in this paper the strength of the 
forcing is adjusted such that the flow remains clearly subsonic, 
that is, the mean Mach number remains below 0.2. Hence for 
all practical purposes the flow can be considered incompressible 
(Dobler et al. 2003). In these simulations no mean flow develops, 
that is t/ = 0, so t/ = H. 

2.4. Simulations 

The relevant equations are solved in a computational domain of 
size L X L X L using periodic boundary conditions. In the case 
of the Roberts flow ( |26] | we fix L by L = 2n/k^). Only four of the 
test field equations (l22l i (those with p - 1, q - c and q = s) are 
solved numerically. With turbulence in the kinematic regime the 
four equations ( |27] i and ( |28] l for U and In p are solved together 
with four of the test field equations (l22T i (again with p - I, q - c 
and q - s). 

Due to the finiteness of the domain in z direction and 
the periodic boundary conditions, quantities like £ and B 
have to be considered as functions that are periodic in z. 
The Fourier integrals used for representing these quantities, 
Q{z) = f Q{k) exp(ife) d{k/2n), turn into Fourier series, Q(z) - 
2 Qn exp(ik„z)/L, where kn = Inn/L and the summation is over 

n = 0, +1, ±2, By this reason only discrete values of k, that 

is k = k„, are admissible in (flJTl-lfTSll. In this framework we may 
determine the a and fjt only for these k„. 

As explained above, the test field procedure yields a and 
fit not as functions of k alone. They may also show some de- 
pendence on z and t. After having averaged over z, time aver- 
ages are then taken over a suitable stretch of the full time series 
where these averages are approximately steady. We use the time 
series further to calculate error bars as the maximum departure 
between these averages and the averages obtained from one of 
three equally long subsections of the full time series. 

In all cases the simulations have been carried out using the 
Pencil Cod43 which is a high-order finite-diff'erence code (sixth 



order in space and third order in time) for solving the compress- 
ible hydromagnetic equations together with the test field equa- 
tions. In the case of the Roberts flow, of course, only the test field 
equations are being solved. 



3. Results 

3.1. Roberts flow 

Let us first recall some findings of earlier work, which are pre- 
sented, e.g., by Radler (2002a,b). We use here the definitions 



(29) 



Adapting the results of analytic calculations in the framework 
of SOCA to the assumptions and notations of this paper (see 
Appendix iBl) we have 



Qf/ffO = '7t/?7tO = 



(30) 



Moreover, in the general case, also beyond SOCA, it was found 
that 



a - aoRtji tp( y/lRm) 



(31) 



with a function (p satisfying 0(0) = 1 and vanishing with growing 
argument. This function has been calculated numerically and is 
plotted, e.g., in Radler et al. (2002a,b). 

Fig. [T] shows results for a and rjt obtained both by general 
test field calculations using (|22] | and under the restriction to 
SOCA using ( |23] l. These results for a agree completely with 
both (l30l l and ( ISTT i. and those for % agree completely with 
(|30] |. Unfortunately we have no analytical results for //t beyond 
SOCA. 

Proceeding now to the a(k) and ^t(^) we first note that in 
SOCA, as shown in AppendixICl 



aoR„ 



1 + (k/kfy 



1 + (k/kf)^ 



(32) 



The corresponding and r]t((), again in SOCA, read 



a(0 



jaokfR^ exp(-fcf 1^1) , 



i/7tofcf^mexp(-fcf|f|).(33) 



[http : //www ■ nordi ta . org/sof tware/penci 1 - code | 



In Fig. 12] results of test field calculations for the functions 
a(k) and fj.ik) with R^ = 10/ V2 * 7.1 are shown. We note 
that rji becomes negative for small k. The same has been ob- 
served with another but similar flow of Roberts type (Radler & 
Brandenburg 2003). For comparison, SOCA results obtained in 
two different ways are also shown: those according to the ana- 
lytic relations ( l32b and those calculated numerically by the test 
field method with ( |23] |. Both agree very well with each other. 

In order to obtain the results for the kernels a(() and f)i(0 
we have calculated integrals as in ( fT4l ) numerically using the data 
plotted in Fig.|2l The results are represented in Fig. [3] Again, an- 
alytical and numerical SOCA results are shown for comparison. 
Note that the profiles of a{0 and f}t{0 beyond SOCA are rather 
narrow compared with those under SOCA, and that of fjtiO sven 
more narrow than that of Q'(^). 

3.2. Isotropic turbulence 

Results for homogeneous isotropic turbulence have been ob- 
tained by solving the hydrodynamic equations ( |27] | and (l28T l si- 
multaneously with the test field equations (l22l l in a domain of 



A. Brandenburg et al.: Scale dependence of alpha effect 



5 



10.0 r 





10 



0.0 0.2 0.4 0.6 



0.8 



1.0 





10.00 






















...■<>■■■■" 




o 
















1.00 














0.10 














0.01 













10 



Fig. 1. Dependences of the normalized a and % on for the 
Roberts flow in the general case, i.e. independent of SOCA (solid 
lines), and in SOCA (dotted lines). 
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Fig. 2. Dependences of the normalized a and fit on k/k{ for the 
Roberts flow with = 10/ V2 » 7.1 (solid lines), compared 
with normalized SOCA results for a/R^ and fi^jR^, which are 
independent of R^ (dotted lines). 



size I? . The forcing wavenumbers kf are fixed hy kfjki -5 and 
10. Instead of the definitions ( |29] ) we use now 



?7tO - \umislkf , Rm - Urms/rjkf . 



Within this framework the dependence of a and rjt on R^ 
has been studied by Sur et al. (2008). They considered two 
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Fig. 3. NormaUzed integral kernels a and % versus kf^ for the 
Roberts flow with R^n - 10/ V2 x 7.1 (solid lines), compared 
with normalized SOCA results for a/Rm and fjt/Rm, which are 
independent of R^ (dotted lines). The full width half maximum 
values of kf^ for a and fji are about 0.5 and 0.2, respectively. 

cases, one with v/rj = 0.1 and another one with Mrms/^^f - 2.2. 
Remarkably, they found that a/ao and 77/770 approach unity for 
^m» 1. 

Figure|4] shows results for a(k) and 77t(^) with v/77 = 1. Both 
a and 771 decrease monotonously with increasing The two val- 
ues of a for a given k/kf but different kf/k[ and R,^ are always 
very close together The functions a(k) and f/tik) are well repre- 
sented by Lorentzian fits of the form 
ao . , , , 77,0 



(35) 



1 + {kiksf ' ' 1-1- {kllkff ■ 

In Fig. |5]the kernels a{() and ^t(^)^ again with v/77 - 1' 
obtained by calculating numerically integrals as in (fT4l i. are de- 
picted. Also shown are the Fourier transforms of the Lorentzian 
fits. 



&(0 



\aoki exp(-^f 1^1) , 



flU) = exp(-2^f 1^1) . 



(36) 



Evidently, the profile of fjt is half as wide as that of a. This cor- 
responds qualitatively to our observation with the Roberts flow 
beyond SOCA, see the crosses in Fig. [3] There is however no 
counterpart to the negative values of 7)1 that occur in the example 
of the Roberts flow. 

The results presented in Figs |4] and |5] show no noticeable 
dependences on R^. Although we have not performed any sys- 
tematic survey in R^, we interpret this as an extension of the 
above-mentioned results of Sur et al. (2008) for a and 771 to the 
integral kernels a and fj^. Of course, this should to be checked 
also with larger values of R^. Particularly interesting would be a 
confirmation of different widths of the profiles of a and 771. 



(34) 4. Discussion 



Our results are important for calculating mean-field dynamo 
models. The mean-field induction equation governing B, here 
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malized wavenumber k/kf for isotropic turbulence forced at 
wavenumbers kf/ki - 5 with R,^ - 10 (squares) and kf/ki = 10 
with Rm = 3.5 (triangles), all with v/rj - 1. The solid lines give 
the Lorentzian fits ( l35T l. 
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Fig. 5. Normalized integral kernels a and fjt versus kf( for 
isotropic turbulence forced at wavenumbers kf/ki - 5 with 
Rm = 10 (squares) and kf/ki = 10 with R^n = 3.5 (triangles), 
all with v/i] - 1. The solid lines are defined by ( [36] l. 

defined as average over x and y, with £ according to dSj, allows 
solutions of the form ^B^) exp(iA:z + /lf)j, Bq, - 0, with the 
growth rate 



A dynamo occurs if A is non-negative. Since a < in all ex- 
amples considered this occurs with the lower sign, and we focus 
attention on this case only. In the limit of a local connection 
between £ and B the ^t(^) and a(k) turn into ^t(O) and a{Q), re- 
spectively. 

When using the definitions ( l29b for the Roberts flow, or (l34b 
for isotropic turbulence, we may write ( l37b in the form 



^ = r]tokf < - 



k a{k/kf)} k 
kf ao J kf 



(38) 



[i] + rjiik)] k^ + a(k)k . 



(37) 



where y - 2 for the Roberts case and y = 3 for the isotropic 
case. Since 1)1 and a depend only via kjkf on k we have chosen 
the arguments kjkf. 

Consider first the Roberts flow, that is dJSl l with y - 2. 
Clearly A is non-negative in some interval < kjki < /cq and it 
takes there a maximum. Dynamos with k/kf > kq are impossible. 
Of course, kq depends on With the analytic SOCA results 
( [32] i we find that /<•() = for small and that it grows mono- 
tonically with R^ and approaches unity in the limit of large R^. 
For small R^ (to which the applicability of SOCA is restricted) 
a dynamo can work only with small k/kf, that is, with scales 
of the mean magnetic field that are much smaller than the size 
of a flow cell. Furthermore, kq never exceeds the corresponding 
values for vanishing nonlocal effect, which is jR^/(l + 5/?^). 
In that sense the nonlocal effect favors smaller k, that is, larger 
scales of the mean magnetic field. With the numerical results be- 
yond SOCA represented in Fig.|2l with R^^ - 10/ VS, we have 
Ko ~ 0.90. ..0.95, again a value smaller than unity. In this case, 
too, a dynamo does not work with scales of the mean magnetic 
field smaller than that of a flow cell. There is no crucial impact 
of the negative values of fj^ for k/kf < 0.8 on the dynamo. 

Proceed now to isotropic turbulence and consider ( |38] | with 
7 = 3. Again A is non-zero in the interval < k/kf < kq and it 
takes there a maximum. Some more details are shown in Fig.|6] 
With the Lorentzian fits i35[ of the results depicted in Fig. |4]we 
find Ko ~ 0.60 for R^ = 10 (and 0.45 for R^ = 3.5). In the 
limit of vanishing nonlocal effects it turns out that kq » 0.82 
for R,^ = 10 (and 0.59 for = 3.5). We have to conclude 
that dynamos are only possible if the scale of the mean magnetic 
field clearly exceeds the outer scale of the turbulence. In addition 
we see again that the nonlocal effect favors smaller k, or larger 
scales of the mean magnetic field. 

These findings may become an important issue especially for 
nonlinear dynamos or for dynamos with boundaries. Examples 
of the last kind were studied, e.g., by Brandenburg & Sokoloff 
(2002) and Brandenburg & Kapyla (2007). In these cases how- 
ever the underlying turbulence is no longer homogeneous and 
therefore the kernels a and fjt are no longer invariant under trans- 
lations, that is, depend not only on ^ but also on z- The finite 
widths of the a and f)t kernels may be particularly important if 
there is also shear, because then there can be a travelling dy- 
namo wave that may also show strong gradients in the nonlinear 
regime (Stix 1972, Brandenburg et al. 2001). 

For another illustration of the significance of a finite width 
of the kernels a and fjt we consider a one-dimensional nonlin- 
ear mean-field model with periodic boundary conditions. We 
modify here the model of Brandenburg et al. (2001, Sect. 6), 
with a dynamo number of 10 (corresponding to 5 times super- 
critical) and R^ - 25, by introducing the integral kernels (l36T l. 
Figure [T] shows the components of the mean magnetic field for 
two different values of kf/ki and for the conventional case where 
the kernels are delta-functions. Note that ki corresponds to the 
largest scale of the magnetic field compatible with the bound- 
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ary condition. It turns out that the magnetic field profiles are not 
drastically altered by the nonlocal effect. Small values of kflk\, 
however, correspond to smoother profiles. 

Let us start again from £ in the form (O, specify there, in 
view of isotropic turbulence, a and according to ( [36] l, and rep- 
resent B(z - and J{z - ^) by Taylor series with respect to A 
straightforward evaluation of the integrals provides us then with 



n>0 



aQ d'"B(z) 77to d'"J(z) 
dz^" {Ikff" dz^" 



)■ 



(39) 



This corresponds to relations of the type (HJ or (O, simply gener- 
alized by taking into account higher than first-order derivatives 
ofB. 

The terms with derivatives of / in ( [39] l can be interpreted in 
the sense of hyperdiffusion. While all of them have the same 
signs in real space, the signs of the corresponding terms in 
Fourier space alternate, which implies that every second term 
acts in an anti-diffusive manner Thus, a truncation of the expan- 
sion should only be done such that the last remaining term has 
an even n, as otherwise ^nti-diffusion would dominate on small 
length scales and cause B to grow beyond any bound. 

There are several investigations in various fields in which 
hyperdiffusion has been considered. In the purely hydrody- 
namic context, Riidiger (1982) derived a hyperviscosity term and 
showed that this improves the representation of the mean veloc- 
ity profile in turbulent channel flows. In the context of passive 
scalar diffusion, Miesch et al. (2000) determined the hyperdif- 
fusion coefficients for turbulent convection and found that they 
scale with n like in Eq. (|39l l. We are however not aware of earlier 
studies differentiating between diffusive and anti-diffusive terms. 

We have investigated the nonlocal cases presented in Fig.|7] 
using truncations of the expansion ( [39] l. However, two problems 
emerged. Firstly, terms with higher derivatives produce Gibbs 
phenomena, i.e. wiggles in B, so the results in Fig.|7]are not well 
reproduced. Secondly, high-order hyperdiffusion terms tend to 
give severe constraints on the maximum admissible time step, 
making this approach computationally less attractive. It appears 
therefore that a direct evaluation of the convolution terms is most 
effective. 



5. Conclusions 

The test field procedure turned out to be a robust method for 
determining turbulent transport coefficients (see Brandenburg 
2005, Sur et al. 2008 and Brandenburg et al. 2008). The present 
paper shows that this also applies to the Fourier transforms of the 
integral kernels which occur in the nonlocal connection between 
mean electromotive force and mean magnetic field, in other 
words to the more general scale-dependent version of those 
transport coefficients. For isotropic turbulence the kernels a and 
fji have a dominant large-scale part and decline monotonously 
with increasing wavenumbers. This is consistent with earlier 
findings (cf. Brandenburg & Sokoloff 2002), where however 
the functional form of the decline remained rather uncertain. 
Our present results suggest exponential kernels, corresponding 
to Lorentzian profiles in wavenumber space. The kernel for the 
turbulent magnetic diffusivity is about half as wide as that for the 
alpha effect. This result is somewhat unexpected and would be 
worthwhile to confirm before applying it to more refined mean 
field models. On the other hand, the effects of nonlocality be- 
come really important only when the scale of the magnetic field 
variations is comparable or smaller than the outer scale of the 
turbulence. 

One of the areas where future research of nonlocal turbulent 
transport coefficients is warranted is thermal convection. Here 
the vertical length scale of the turbulent plumes is often com- 
parable to the vertical extent of the domain. Earlier studies by 
Miesch et al. (2000) on turbulent thermal convection confirmed 
that the transport of passive scalars is nonlocal, but it is also more 
advective than diffusive. It may therefore be important to also 
allow for nonlocality in time. This would make the expansion 
of passive scalar perturbations more wave-like, as was show by 
Brandenburg et al. (2004) using forced turbulence simulations. 
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Appendix A: Justification of equations (H) and ([12) 

In view of Q we start with equation (O for £>, put £o - 0, 
assume that K o B is a purely spatial convolution. Applying 
then the Fourier transform as defined by Q(x) = j Q{k) exp(i^ ■ 
x) d^ik/ln) we obtain 



8i(k) = Kij(k)Bj(k) . 



(A.l) 



Since fi and B have to be real we conclude that K*j{k) = Kij(-k). 

Further the assumption of isotropic turbulence requires that Kij 
is an isotropic tensor. We write therefore 

Kij = aik)Sij + a{k)kikj + \i)i{k)eijkkk (A.2) 
with a, a' and i), being real functions of = |^|. Considering 
further that k ■ B - and ik x B - J v/e find 



Gik) = aik)B{k) - iit(k)J{k) . 



(A.3) 



Transforming this in the physical space we obtain immediately 
©. 

In view of (fT2] i we start again from equation (|3]l and put £>o - 
but we have to consider KoB now as a convolution with respect 
to z only. Applying a Fourier transformation defined by Q{z) - 
j Q{k)sxpiikz)d(k/27T) we obtain a relation analogous to (lA.ll) . 



Si(k) = Kij(k)Bj(k) . 



(A.4) 



and may now conclude that K*j{k) - Kij{-k). We arrive so at 
Kij = aij(k) + ikfi'ijik) (A.5) 
with real tensors and fj'.j, which are even in k. Combining 
jA.Aj and jA.Si and considering that the ikBi can be expressed 

by the 7, {ikBi - J2, ikBi = -Ji, ikB^ - 0) we may confirm 
first ( fT3] ) and so also ( fT2b . 

Appendix B: IVIean-field results for the Roberts flow 

A mean-field theory of the Roberts dynamo, developed in view 
of the Karlsruhe dynamo experiment, has been presented, e.g., in 
papers by Radler et al. (2002a,b), in the following referred to as 
R02a and R02b. There a fluid flow like that given by ( |26] | is con- 
sidered but without coupling of its magnitudes in the jcy-plane 
and in the z-direction. The mean fields are defined by averaging 
over finite areas in the xy-plane so that they may still depend on 
X and y in addition to z. As shown in the mentioned papers £, 
when contributions with higher than first-order derivatives of B 
are ignored, has then the form 

£ = -a^[B-{z-B)z]-Pi.'VxB-%-j3j[z-{VxB)]z 

-/3iZx[WzB) + (z-V)B] (B.l) 

with constant coefficients aj_, /3±, j3\\ and ^63 [see (R02a 9) or 
(R02b 9)]. Reducing this to the case considered above, in which 
B depends no longer on x and y, we find 



where V xB - zx dB/dz, and 

a = Tli^P^+P^- (B.3) 

Results for a^, /3j_, /3\\ and fij, obtained in the second-order 
correlation approximation are given in (R02a 19) and (R02a 49) 
as well as in (R02b 19) and (R02b 38). When fitting them with 
Uj_ — (2/n)uo, M|| = V2(2/7r)^Mo, ^r/a — ki, V^tt/o = kf, R^j^ - 
V2/?m and /?,n|| = (8/7r)7?n, to the assumptions and notations used 
above we find just Likewise (R02a 20) and also (R02b 20) 
lead to (I3B- 

Appendix C: a and 771 under SOCA for Roberts flow 

Let us start with the relation jB.lj and subject it to a Fourier 
transformation so that 



(C.l) 



£ = a[B - (z ■ B)z] - 77tV X B , 



(B.2) 



£ = mxA = q'[B-(z- B)z] -ikfjtZxB. 
From the induction equation we have 

77(V2 - k^)b = - (V H- i;tz) x{uxB), kk^O. (C.2) 
The solution b reads 

b = -^5^{zxV('^- ViA)-Af(^-VtA)z 

+ik[z xV^(z-% + knlj{t - (z ■ ^)z)]} , (C.3) 

where u is according to ( |24] | expressed by ijj. Calculate now £ v or 
&y and note that 1]/^ - jimlky)^ and (dift/dx)^ = - if/d^ift/dx^ = 
4m?. When comparing the result with ( IC.lb we find immediately 
Using then relations of the type (fT4l) we find also ( l33T l. 
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